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J> ■ ABSTRACT 

QQ ^ We present 70 and 160 /dm Herschel science demonstration images of a field in the Orion A molecular cloud that contains the 

prototypical Herbig-Haro objects HH 1 and 2, obtained with the Photodetector Array Camera and Spectrometer (PACS). These 

^N^ ■ observations demonstrate HerscheFs unprecedented ability to study the rich population of protostars in the Orion molecular clouds at 

• ■ the wavelengths where they emit most of their luminosity The four protostars previously identified by Spitzer 3.6-40 fim imaging 

1^ I and spectroscopy are detected in the 70 yt/m band, and three are clearly detected at 160 yt/m. We measure photometry of the protostars 

^? , in the PACS bands and assemble their spectral energy distributions (SEDs) from 1 to 870 fim with these data, Spitzer spectra and 

^p ■ photometry, 2MASS data, and APEX sub-mm data. The SEDs are fit to models generated with radiative transfer codes. From these 

fits we can constrain the fundamental properties of the protostars. We find luminosities in the range 12 - 84 L© and envelope densities 

spanning over two orders of magnitude. This implies that the four protostars have a wide range of envelope infall rates and evolutionary 

states: two have dense, infalling envelopes, while the other two have only residual envelopes. We also show the highly irregular and 

S^ ■ filamentary structure of the cold dust and gas surrounding the protostars as traced at 160 yt/m. 

^^ ' 

^ ■ Key words. Stars: formation — Stars: protostars — circumstellar matter — Infrared: ISM — Infrared: stars 

1, Introduction peak of the spectral energy distributions (SEDs) of protostars, 

, ^. t t t , , . ' n which are dominated by thermal emission from a cold (-10 K) 

The Orion molecular clouds are the most active region of star e^^eiope Measuring the peak of the SED allows firm estimates 

formation within 500 pc of the Sun, where the Spitzer Space ^f j^^ bolometric luminosities and envelope densities of the pro- 

Telescope identified over 400 likely protostars in the Orion A and tostellar svstems 

B clouds (Megeath et al., in prep). The region is home to both Wit h the Photodete ctor Array Camera and Spectrometer 

clustered and distributed star formation and hosts both high- and (^^g; iPoditsch et alJ I20T0I) aboard Herschel dPilbratt et alJ 

ow-mass protostars The Herschel Space Observatory s capa- |^^ ^^ j^^^^ obtained 70 and 160 m images of a field in 

bihties in the far infrared are crucial for sampling the expected ^^^ Lynds 1641 region. The field contains the intermediate-mass 

— — — , .... Herbig B9e star V380 Ori ( Hillenbrand et alJ il992). 28 infrared 

Herschel is an ESA space observatory with science instruments -j ^-c ju u »■ vu o -^ ^n/r tu 

., ,, „ , jfl- • IT .• . .• J -.1 • excess sources identified by observations with 5»ifzer (Megeath 

provided by European-led Principal Investigator consortia and with im- , . , , r n r, . • > . • 

portant participation from NASA. This work includes data acquired ^^ al, m prep.), and a variety of outflow phenomena mcludmg 

with the Atacama Pathfinder Experiment (APEX; E-082.F-9807, E- the well-known HH 1 and 2 (e.g., Bally et alJ|2002D and > 8 

284.C-5015). APEX is a collaboration between the Max-Planck-Institut protostellar outflows (e.g., Stanke et al. 2002). This is the sci- 

fur Radioastronomie, the European Southern Observatory, and the ence demonstration field for the //^r^c/z^/ Orion Protostar Survey 

Onsala Space Observatory. (HOPS), a 200-hour open-time key program that will obtain 
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PACS imaging of 133 fields, 5' to 8' in diameter, containing 278 
protostars and PACS spectroscopy of a subset of 37 protostars. 

Here, we present photometry of the four protostars in the 
Herschel field at 70 and 160 yum and combine these data with 
S pitzer and ground-base d data. With the radiative transfer code 
of IWhitney et al.l (l2003b . we generate model SEDs and find that 
the four protostars exhibit a large range of luminosities (12 < 
LJLq < 84) and envelope densities spanning over two orders of 
magnitude. This implies that two protostars have dense, infalling 
envelopes, while the other two have only residual envelopes. 



2. Observations and Data Reduction 

An 8' square field with central coordinates a = 5^36'"22^.05, 
S = -6M5'4r'.23 (J2000) was observed on 2009 October 9 (ob- 
serving day 148; observation IDs 1342185551 and 1342185552) 
in the 70 jim ("blue") and 160 jim ("red") bands available with 
PACS, which have angular resolutions of 5.2'' and 12'', respec- 
tively. We observed our target field with homogeneous coverage 
using two orthogonal scanning directions and a scan speed of 
20"/s. Each scan was repeated 5 times for a total observation 
time of 1468 s per scan direction. The eff'ective sampling rate of 
the detectors is 10 Hz. The data were processed from raw teleme- 
try to final images with the Herschel Common Software System 
(HCSS) version 3.0 build 919, using version 4 of the flux calibra- 
tion files Q We fol lowed the standard pro cessing steps for PACS 
data described by lPoglitsch et al.l (12010 ) with these exceptions: 
We identified cosmic rays for each spatial sky pixel as those val- 
ues which were larger than 10 standard deviations from the mean 
signal. Several extraneous calibration measurements were inter- 
spersed with the HOPS target observations. These were masked 
and removed from the data cube with an additional 430 readouts 
following each of these calibration measurements to mask signal 
drifts induced by the calibration source. 

After the initial processing, the two orthogonal scan observa- 
tions were combined for the final mapmaking step. Two diflTerent 
mapping approaches are used for this purpose: Method 1 is used 
exclusively for point source photometry, while Method 2 is used 
to display images. 

Method 1: Mapping with local sky subtraction: First, we re- 
move the signal drifts (whether correlated or due to the 1// 
noise) by subtracting a local "sky" value from each readout from 
each bolometer pixel. The local sky is estimated as the median 
value within a window of size ±20 readouts. The final mosaic 
is then created by spatially averaging all overlapping bolome- 
ter pixels using the HCSS routine "photProject." To protect the 
integrity of the point source PSF, all readouts within 20" of a 
point source are ignored during the sky median calculation. This 
processing preserves all point and compact sources in the image 
and provides the proper photometry comparison between HOPS 
target objects and the flux calibration standards, which use the 
same reduction scheme. However, this processing removes all 
emission at spatial scales larger than the median window size. 

Method 2: Mapping without local sky subtraction: Method 1 
is necessary only for accurate photometry of point sources. We 
also create maps by removing only the pixel-to-pixel electronic 
oflTsets in PACS images, using the median value of the entire 
time stream of a single pixel to estimate its oflfset signal value. 
Unlike Method 1 , this approach does not remove the (spatially) 
extended emission. However, it also does not mitigate the 1// 




^ HCSS is a joint development by the Herschel Science Ground 
Segment Consortium, consisting of ESA, the NASA Herschel Science 
Center, and the HIFI, PACS, and SPIRE consortia. 



Fig. 1. Three-color composite image of the HH 1-2 region. Blue is 
NEWFIRM K,, green is PACS 70 yum, and red is PACS 160 fim. 



drifts, which add so-called "striping" or "banding" in the final 
maps. As for Method 1, we use the "photProject" HCSS routine 
to spatially coadd individual array readouts for mapmaking. 

A Ks image of the field was acquired with NEWFIRM, the 
NOAO Extremely Wide Field Infrared Imager, on the KPNO 
4 m telescope, and the data were redu ced with the NOAO 
NEWFIRM Pipeline (IS waters et al.ll2009l) . The on-source time 
was 1 1 minutes over most of the field of view. Images at 350 and 
870 Jim were acquired at APEX with SABOCA and LABOCA, 
respectively. The observing and data reduction proce dures for 
the APEX images are described in lStanke et al.l (l2010l) . 



3. Results 

3.1. Imaging and Photometry 

Figure [T] shows a composite of the final map created using 
Method 2 for the 70 and 160 yum PACS channels. Figures [2] and 
[3l available in electronic form only, show the separate 70 and 
160 yum images and are annotated with source names. 

The bright blue source in the north of the field is the reflec- 
tion nebula NGC 1999; the dark tri-lobed fea ture seen toward 
this nebula is discussed in Stanke et al. (2010). In the center of 
the image is a triangular arrangement of protostars. Here we 
use their designations for the HOPS program: 165, 166, 168, 
and 203. HOPS 166 (HH 147 MMS; Chini et al. 2001) is the 
relatively isolated source at the northeastern corner of the tri- 
angle, HOPS 168 (HH 1-2 MMS 2) is at the western corner, 
and HOPS 165 and 203 are the pair of overlapping sources 
(separated by 13'') at the southern corner. HOPS 203 (HH 1- 
2 MMS 1), the brighter of the pair in the PACS bands, is the 
source of the HH 1-2 outflow and contains the radio sources 
VLA 1 and 2 (Rodrig uez et al.|[l99Q|) . IChini et al.i report an ad- 
ditional source HH 1-2 MMS 3, 22'' southwest of HOPS 168, 
that corresponds to extended emission at 160 yum with no ap- 
parent point source at 70 yum. Falling nearly along the line be- 
twee n HOPS 168 and 203 is the C-S star, a classical T Tauri 
star (ICohen & Schwartzlll979h . To the southeast of HOPS 203 at 
a = 5^36^25^.3, S = -6°47'18" is a knot of emission presum- 
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ably shock heated by the HH 2 outflow. At 160 yum, only HOPS 
166, 168, and 203 appear, while HOPS 165 is not detected. The 
160 yum band also traces cold dust in the surrounding cloud ma- 
terial, showing an irregular, filamentary structure. 

PACS photometry of the four protostars appears in columns 
9 and 10 of Table [1] We obtained simple aperture photometry for 
the relatively isolated protostars: HOPS 166 and 168 in the blue 
and red bands and HOPS 203 in the red band. In these cases, 
we used a 16'' aperture with subtraction of the median signal in 
a background annulus extending from 18'' to 22". The results 
were corrected according to measurements of the encircled en- 
ergy fraction provided by the PACS consortium (priv. comm.). 

For the HOPS 165/203 pair at 70 yum, point- spread function 
(PSF) fitting was required to separate the fluxes of the two proto- 
stars. We fit the fainter HOPS 165 with a PSF constructed from 
observations of Vesta (PACS consortium, priv. comm.). Aperture 
photometry for HOPS 165 was performed on the best-fit PSF, 
and aperture photometry for HOPS 203 was performed on the 
data after subtraction of the HOPS 165 model. At 160 yum, we 
report an upper limit for HOPS 165; this is the largest flux den- 
sity for which a model PSF can be added at the source position 
before it appears as an asymmetry in th e HOPS 203 image. 

According to Poglitsch et al.l (l2010h . the calibration accuracy 
for PACS is within 10% in the blue band and better than 20% in 
the red. The formal uncertainties associated with each source 
(i.e., the RMS of the signal in the sky annulus) are much less, 
< 1%, except for the case of HOPS 165, where fitting a point- 
spread function to a faint source yields a 10% uncertainty. 

The PACS photometry data are supplemented by Spitzer 
IRAC and MIPS photometry (Megeath et al., in prep.), Spitzer 
IRS spectros copy, and APEX SA BOCA and LABOCA sub-mm 
photometry (IStanke et al.ll2010h . For HOPS 166, near-infrared 
J/H/K photometry was available from the Two Micron All Sky 
SurveyQ The Spitzer positions and 3.6 - 870 yum photometry for 
the HOPS protostars appear in Table [T] Systematic uncertainties 
are given in a note to the table. 



3.2. SED Modeling 

We u se a Monte Carlo radiative transfer code (IWhitney et al.l 
l2003h to calculate model SEDs for the four protostars. The code 
features a central star and flared disk, which emit photons that 
can then be scattered or absorbed and re-emitted by dust in ei- 
ther the disk or an envelope. The envelope den s ity is defined 
by the rotating collapse solution of iTerebey et al.l (Il984l) . plus a 
bipolar, evacuated cavity. 

We use the same dust model as iTobin et aP (l2008l) . which 
contains larger dust grains than a standard ISM dust model. The 
grain size distribution is defined by a power law n (a) oc a~^-^, 
with 0.005 yum < a < I yum. We use dust grains composed of 
graphite ^graph = 0.0025, silicates ^sii = 0.004, and water ice 
fice = 0.0005; abundances (^) are relative to gas and imply a 
gas to dust ratio of 133. Our sub-mm opacities exceed those 
of the well-known Mi l ky Wa y Case B (Ry = 5.5) mixture of 
IWeingartner & Draina (l200lh by a factor that reaches a maxi- 
mum of 5 at 600 yum. 

The model parameters are set to typical values for low-mass 
protostars; we fit the SEDs by varying seven of the m: the system 
luminosity L, the reference envelope density pi (iKenyon et al.l 




Fig. 4. Observed SEDs (black) with best-fit models (blue). Plus signs 
indicate 2MASS, Spitzer, and APEX photometry; the open triangle is an 
upper limit. Black curves show Spitzer spectra. Filled symbols indicate 
Herschel photometry: diamonds represent detections, while the filled 
triangle is an upper limit. 



^ The Two Micron All Sky Survey (2MASS) is a joint project of the 
University of Massachusetts and the Infrared Processing and Analysis 
Center/California Institute of Technology, funded by NASA and the 
National Science Foundation. 



Il993l) . the outer radius of the envelope 7?env, the opening angle 
of the envelope cavity ^cav, the mass of the disk Mdisk, the in- 
clination angle /, and the foreground ex tinction Ay. (Fo reground 
extinction is applied with the laws of lMcClur3 l2009l suitable 
for star- forming regions.) In fitting the sources, we emphasize 
the mid to far IR over the near IR, since the near IR is highly 
dependent on the scattering properties of the dust, the geometry 
of the inner disk, and the geometry of the outflow cavity. Thus, 
when fitting a source we first adjust the luminosity and density 
to get the best fit to the mid to far IR, then we find the best com- 
bination of cavity opening angle, inclination, and (if necessary) 
foreground reddening to fit the 10 yum absorption feature and the 
near-IR emission. In general, the fits are insensitive to 7?env and 
Mdisk- However, for HOPS 165, it was necessary to adjust these 
two parameters. The best-fit parameters were determined by vi- 
sual comparison of the models to the observed SEDs and are 
listed in Table |2l 

The models are compared to the photometry and spectra in 
Figure m For these models, 4 x 10^ photons were run through 
the Monte Carlo code. The code generates output for apertures 
ranging from V to 16'' in one-arcsecond steps, and the choice 
of aperture for the plotted SED varies with wavelength, as given 
in the note to Table [T] An interpolation scheme bridges the gaps 
between disparate apertures. We assume a distance of 420 pc 
dMenten et al.ll20Q7i) . 

4. Discussion 

The Herschel observations have detected four protostars in the 
HH 1-2 region. All were previously identified with Spitzer 
(Megeath et al., in prep), but Herschel has cleanly separated 
the sources for the first time in the far IR, allowing accurate 
photometry. From the luminosities and densities in Table [31 we 
estimate the infall rate from the envelope onto the central star- 
disk system, the luminosity due to accret ion onto the star, and 
the evolutionary state for each protostar. iKenyon et al.l (Il993h 
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Table 1. Protostar Photometry 



HOPS 


RA (J2000) 


Dec. (J2000) 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


[24] 


[70] 


[160] 


[350] 


[870] 


Source 


(hms) 


r ' ") 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


165 .. 


5 36 23.54 


-6 46 14.6 


0.014 


0.052 


0.11 


0.13 


0.64 


1.1 


<0.3 


<3 




166.. 


5 36 25.13 


-6 44 41.8 


0.66 


0.83 


0.97 


1.17 


4.65 


10.9 


11.1 


4.6 


0.33 


168.. 


5 36 18.93 


-6 45 22.7 


0.0077 


0.030 


0.041 


0.038 


3.66 


87.3 


87.7 


24 


0.94 


203.. 


5 36 22.84 


-6 46 06.2 








0.0091 


0.54 


26.6 


75.7 


28 


1.2 



Notes. Aperture radii are 3'' from 1.2 to 2.2 yum (not shown), 2" from 3.6 to 8.0 //m, 5'' at 24 yum, and 16" at 70 yum and longward. Systematic 
uncertainties are 3% from 1.2 to 2.2 yum, 5% from 3.6 to 24 yum, 10% at 70 yum, 20% at 160 yum, 40% at 350 yum, and 20% at 870 /um. 



Table 2. Adopted ]V[odel Parameters 



HOPS 


L 


Pi 


-f^env 


^cav 


Mdisk 


/ 


Av 


Source 


(Lq) 


(g cm 3) 


(AU) 


n 


(Mo) 


n 




165.. 


12 


7.5 X 10-^' 


10^ 


30 


1 X 10-^ 


20 


35 


166.. 


23 


1.5x10-15 


10^ 


25 


5 X 10-2 


40 


4 


168 .. 


84 


3.0 X 10-13 


10* 


40 


5 X 10-2 


75 





203.. 


23 


2.6 X 10-13 


10* 


40 


5 X 10-2 


75 






The envelope density at 1 AU in the limit of no rotation. 



show that for a protostellar envelope, the infall rate is Menv = 
1.9 X 10"^ (pi/10-^5 g cm"^) {MJUq^I^ Mq yx-\ The accre- 
tion luminosity can then be written as Lace = GM^Mdisk/R* = 
5.9 (pi/10-15 g cm-3) (M./Mo)3/2 (RJRq)-' (Md,sk/Menv) Lq, 
where Mdisk is the accretion rate from the disk onto the star. 
To estimate Mgnv and Lace, we (initially) assume Mdisk = Menv, 
and we adopt ste llar radii, luminosities, and masses from the 
ISiess et aP (l2QQ0h online models of pre-main- sequence stars at 
an age of 3 x 10^ yr. The final conclusions are not sensitive to 
the exact stellar parameters chosen. 

HOPS 166 is modeled as a luminous star-disk system with 
a low-density envelope seen through a few magnitudes of vi- 
sual extinct ion. (The modeled i nclination of 40° is considered a 
lower limit; l Eisloffel et al.lll99 4 find that the outflow associated 
with this source is close to the plane of the sky.) The best-fitting 
central star from the Siess et al. models has a mass of 2.2 Mq, 
implying an envelope infall rate of 4 x 10"^ Mq yr"^ and an ac- 
cretion luminosity that is 20% of the total luminosity. The low 
accretion luminosity and low envelope mass (inferred from the 
model parameters) of only 0.02 Mq imply that HOPS 166 is in 
the late stages of protostellar e volution and that t he central star 
has accreted most of its mass. (IChini et al.l 12001 classified this 
source as a deeply embedded Class object based on SCUBA 
and IRAM mapping at 450, 850, and 1300 yum.) 

In contrast, HOPS 168 is much more embedded and lumi- 
nous than HOPS 166. Its envelope mass is 2.7 Mq. The im- 
plied stellar mass is 0.3 Mq, the envelope infall rate is 3 x 
10"^ Mq yr"\ and the accretion luminosity is more than 95% 
of the total. A star more massive than 0.3 Mq is possible if 
Mdisk < Menv, meaning inf ailing r natter is piling up on the disk, 
leading to episodic accretion (e.g. JVorobyov & Basull2QQ5h . For 
example, the central star could have a mass as high as 1.8 Mq if 

Mdisk =0.1 Menv. 

The two remaining protostars, HOPS 165 and HOPS 203, 
are separated by only 13'', or a projected 5500 AU. The SED 
of HOPS 165 drops off' precipitously beyond 30 fim. This re- 
quires a very small, tenuous envelope and a low-mass disk. The 
ffux from the moderately luminous star-disk system is seen be- 
hind 35 magnitudes of visual extinction. Our interpretation is 
that HOPS 165 must be seen through the dense envelope of the 



nearby HOPS 203 (Menv = 2 .4 Mq). HOPS 203 itself is a ^ 
binary (iRodriguez et al.ll 19901) . If the proximity of HOPS 165 is 
not due to chance, this region is home to a hierarchical multiple 
system of three protostars within a projected radius of 5500 AU. 
Accordingly, the small envelope size of the HOPS 165 model 
may result from its proximity to HOPS 203. The implied stel- 
lar mass of HOPS 165 is 1.4 Mq, the envelope infall rate is 
2 X 10"^ Mq yr"\ and the accretion luminosity is 10% of the 
total. On the other hand, the implied stellar mass of HOPS 203 
is 0.1 Mq, the envelope infall rate is 2 x 10"^ M© yr"\ and the 
accretion luminosity is more than 95% of the total. Again, the 
central star may have a higher mass if Mdisk < Menv We assume 
that the accretion is dominated by one member of the 3'' bi- 
nary, but the results will not change significantly if both accrete 
equally. The 160yum PACS measurement for HOPS 203 exceeds 
the fit by a factor of 2; this may be due to cold envelope material 
in our aperture that is not accounted for in our models. 

We conclude that two of the protostars (HOPS 168 and 203) 
are in an active state of mass infall and accretion, while the other 
two (HOPS 165 and 166) have only residual envelopes. This 
finding demonstrates HerscheVs unique and critical contribution 
to the audit of the ffow of mass from the outer protostellar enve- 
lope onto the central protostar. 
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Fig. 2. PACS 70 nm image of the HH 1-2 region. 
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Fig. 3. PACS 160 yum image of the HH 1-2 region. 



